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We show that the non-local polarization response in a multimode circuit-QED setup, devised from 
a Cooper pair box coupled to a long chain of Josephson junctions, provides an alternative route to 
access the elusive Kondo screening cloud. For moderate circuit impedance, we compute analytically 
the universal lineshape for the decay of the charge susceptibility along the circuit, that relates 
to spatial entanglement between the qubit and its electromagnetic environment. At large circuit 
impedance, we numerically find further spatial correlations that are specific to a true many-body 
state. 


Superconducting circuits constitute at present one of 
the most versatile platforms for quantum engineering, 
due to their macroscopic - yet fully coherent - building 
blocks p], offering great possibilities in tunability and 
design [2], Beyond their promising use in quantum in¬ 
formation processing, such large-scale electrical devices 
made out of Josephson junctions can be viewed as meta¬ 
materials 0 0 where light-matter interactions can be 
explored in uncharted territory Baca . These systems al¬ 
low for unusual physical regimes where the effective fine 
structure constant can become of order one mm , so 
that ultra-strong coupling between a single two-level sys¬ 
tem and a large number of environmental modes leads to 
a vacuum with non-trivial entanglement properties [12] . 
Our aim in this Letter is to develop a simple physical pic¬ 
ture for this massively entangled state, that we dub the 
Josephson-Kondo screening cloud, and to propose a real¬ 
istic setup where its subtle correlations can be unveiled 
experimentally. 

The use here of the Kondo terminology m, usually 
associated to the quenching of a magnetic moment in 
a Fermi sea (another example of a complex many-body 
vacuum) may come as a surprise. However, deep con¬ 
nections between this purely electronic phenomenon and 
strongly dissipative two-level systems have been known 
for decades (LH - fTH] . a property which we will exploit 
here and further elaborate on. In fact, one challenging 
and still open question in the Kondo realm concerns pre¬ 
cisely the evidence for the magnetic screening cloud, or 
in other words, spatial fingerprints of the entanglement 
structure between the localized spin and its surrounding 
electron bath [T0T121 1 . Despite numerous experimental 
efforts mm and innovative theoretical proposals [251 - 
129] . two important difficulties are typically encountered. 
First, measuring weak long-range magnetic correlations 
is a daunting task, and an all-electrical setup would be 
much preferred. Second, spatial Kondo correlations show 
rapid 2 kp Friedel oscillations due scattering processes 
across the Fermi surface, that complicate the signal to 
analyze Ell¬ 
in order to bypass these difficulties, we propose in 



FIG. 1. (color online) Schematic view of the proposed su¬ 
perconducting circuit supporting two dressed qubits states, 
where one or two Cooper pairs in excess on the leftmost island 
are each accompanied by a different charge polarization cloud. 
Josephson tunneling between the two configurations generates 
entanglement between the qubit and its environment, leading 
to a state that possesses the same singlet structure as the 
magnetic Kondo screening cloud, but in the charge, rather 
than the spin sector. 


Fig-Ha simple all-superconducting circuit. Here, one end 
of a long chain of Josephson elements is coupled to a sin¬ 
gle superconducting quantum dot with a large Coulomb 
charging energy. The dot is tuned to a charge degeneracy 
point, so that the state space of the dot is reduced to two 
active charge levels that contain an excess of either one or 
two Cooper pairs. This realizes the analog of a quantized 
spin 1/2 moment, but in the charge sector, which is better 
suited for electrical control with local gates. The low en¬ 
ergy modes of the environment are plasmons that couple 
both capacitively and inductively (via Josephson tunnel¬ 
ing) to the qubit. Quenching the charge qubit in this 
superconducting environment is preferable to quenching 
with a normal electron gas because, unlike the capaci¬ 
tive couplings between elements in the superconducting 
array, electron-electron interactions in a normal electron 
gas act to distort the charge screening cloud. 

The interplay of charging and inductive couplings leads 
to the following physics. The capacitive interaction be- 
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tween the charge qubit and a high impedance bath gener¬ 
ates the dressing of each of the two qubit states by a dis¬ 
tinct charge polarization cloud, as depicted in Fig. [I] The 
Joseplison coupling, rooted in quantum tunneling, forces 
in first approximation a quantum superposition of these 
two dressed qubit states, leading to non-trivial entangle¬ 
ment between the qubit and the bath. The spatial corre¬ 
lations of this polarization cloud can then be probed by 
modulating a comb of local gates applied along the chain, 
while recording the charge response of the two-level sys¬ 
tem. Qubit measurements can be performed by weakly 
coupling the qubit to a superconducting resonator. We 
will show that such measurements allow one to map pre¬ 
cisely the usual magnetic Kondo screening cloud correla¬ 
tions, averaged over fluctuations on the scale of the Fermi 
wavelength. This bypasses the issue of the 2oscilla¬ 
tions mentioned above. A final key aspect pertaining 
to Kondo physics is universality, namely that phenom¬ 
ena beyond a model-specific short distance cut-off are 
insensitive to microscopic details. At low to moderate 
impedance, we will analytically compute the full spatial 
dependence of the Josephson-Kondo screening cloud, and 
demonstrate its universality beyond the first few sites 
of the array. For large dissipation (circuit impedance 
roughly larger than half the resistance quantum), we will 
numerically calculate the universal cloud, using a recently 
developed technique based on a generic coherent state de¬ 
scription of environmental wavefunctions, and find that it 
only depends on the dissipation strength and the Kondo 
length. 

Model. At temperatures well below the superconduct¬ 
ing transition temperature, the proposed Josephson cir¬ 
cuit is fully governed by the dynamics of conjugate phase 
<pi and charge n* degrees of freedom on the various nodes 
i of the array, and obeys to the following Hamiltonian 


the chain is not small. This produces the anharmonicity 
in the second line of Eq. ([!]) , which is responsible for the 
non-trivial physics that we discuss now. 

A deep connection to quantum optics can be made fol¬ 
lowing three standard steps HU, starting with a normal 
mode diagonalization of the quadratic part in Eq. ([l]), fol¬ 
lowed by a unitary transformation H = UHU' , with U = 
e z{n a -no)4> i, anc j fi na iiy a truncation of the Cooper pair 
box Hilbert space (at site * = 0) to two nearly degenerate 
states, owing to the parametrically large charging energy 
at the end of the chain, namely (2 e) 2 /(C gd + C d ) E Jd 
(see Supplementary Material [315] for details). Using the 
replacement with Pauli matrices (no — no) -± <j z /2 and 
e *0o CT + ; anc i sending IV —> oo, this gives readily the 
Hamiltonian: 


H= dk 


Ukb\b k - g k {b\ + 4)yf] - E Jd~ c f, (2) 


2 J 


Wfc = 2sin(fc/2)^ 

1 C gd 


(2e) 2 Ej 


C g + 4C sin 2 (fc/2) ’ 


(3) 


Qk — 


Wk 


y/2 Cgd + Cd sin(A;/2) y 2irEj 


Wfc 


cos[fc/2 - 4], (4) 


the phase shift obeying sin(4) = sin(4 — k), with 
SI = CdC g d/[Cg[C d + Cgd)]- Under the assumptions 
made before on the magnitude of the capacitances, one 
finds that 61 -C 1, so that the phase shift is small in 
practice. For low frequencies w -C -J ( 2e) 2 Ej/C , we find 
a linear spectral density of modes J(u>) = 2irauj , with 
effective fine structure constant (see Supplementary Ma¬ 
terial [39] for a discussion on this interpretation of a) 


1 ( Cgd V 2 /W 
2tt \ C g d + C d J V EjC' 


(5) 


H = ( ni ~ {P Oi ^ 

i,j=0 

E Ar " 1 

+ -p -EjdCOS^o -< 4 ), ( 1 ) 

i=1 

with hi the average number of Cooper pairs on island i. 
N +1 the total number of islands, and the commutation 
relation [<l>j,nj>] = i6jy (we set H = 1 in what follows). 
The matrix C contains on-site capacitances Coo = C gd + 
C d , C\\ = Cg + C d + C, Ca = Cg + 2C', i > 2, on 
the diagonal, as well as nearest neighbor capacitances 
Coi = Cio = —C d , Cii+i = C t+li = — C, i > 1 |32j. All 
other entries of C are zero. We assumed that the charging 
energy in the chain (for sites i > 0) is much smaller 
than the Josephson couplings, namely (2e) 2 /(2C'+C g ) <C 
Ej , so that phase differences in the chain are small. We 
therefore expanded the Josephson terms —Ej cos(<4 — 
<j>i+i), i > 1, to quadratic order. We stress that the phase 
difference between site i = 0 (the dot), and site i = 1 of 


Protocol. The observable that characterizes the spatial 
profile of the screening cloud is 

Xj = {(no ~ no)(nj ~ fij)) - ( 6 ) 

which quantifies the difference in polarization on site 
j of the chain, associated with the distinct charge 
states no — no = ±1/2 of the qubit. This quantity 
can be extracted from frequency-dependent linear re¬ 
sponse measurements as follows. We confine ourselves 
to zero temperature, as a finite temperature simply in¬ 
troduces exponential decay beyond the thermal wave¬ 
length [35]. A small AC gate voltage V(t) = 6Ve lult 
is applied at the site j of the chain, causing a pertur¬ 
bation 6H = 2eC g (C~ 1 )jjSVe lult (nj — fij). From the 
Kubo formula, the linear response of the qubit charge 
at site 0 is (no — ho) = 2eC g (C~ 1 )jjSVx(j, w) with 
X(j, u) = —i / 0 + °° dte lult ([no(t), rij(0)]). At zero temper¬ 
ature if follows from the fluctuation-dissipation relation 
that Xj = - r duJ Im x(j, w)/tt. 
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After completing the same transformations on Xj used 
previously to map Hamiltonian (jTJ) into |2]), one finds 
from (|6| the general expression for the Josephson-Kondo 
cloud: 


i r 

Xj = 2 J Q dk ° jk 




( 7 ) 


0jh = 4 


Ej 

2nu>k 


sin(fc/2) cos[/u(j — 1/2) — Sk\- (8) 


Regime of intermediate fine structure constant. For 
up to moderate values of the circuit impedance (typi¬ 
cally a < 0.2), it was shown recently ]T2lf33lf35] that the 
full ground state |'F) describing the coupled qubit and 
Joseplison chain system assumes with excellent accuracy 
a simple form I’F) = (|/) |f) + |-/) It)) /v^, where 


l±/) = exp 


± / dk f k (b\ — b k ) 
Jo 


| 0 ), ( 9 ) 


and the oscillators are displaced by: 




2 u>k + A R ’ 


A R = E Jd exp 


-2 


dkfl 


■ (10) 


This provides the mathematical foundation for the phys¬ 
ical picture shown in Fig. [l] underlying the singlet-like 
structure of the wavefunction, that leads to entanglement 
between the qubit and environmental degrees of freedom. 

Substituting the ground state into the Josephson- 
Kondo cloud polarizability ([?]), one finds for j > 1: 

_ A R C gd rdk cos(fc/2) cos [k(j - 1/2) - 6 k \ 
C g d + Cd Jo 2tt Wfc + A R 

(n) 

This expression depends on all microscopic parameters 
of the device. However, for j sufficiently larger than 
VCjCg, one finds universal single-parameter scaling: 
(See the supplementary material [39] for further detail.) 


Xj = — Z Re < e B o J r 


n • A-R • 
0 


( 12 ) 


where r(a, z) is the incomplete Gamma function, Eq = 
y/(2e) 2 Ej/C g , and 


Z" 1 = 2tt 1 


C d \ E, 


Cgd) A R 


(13) 


The emergence of a universal scaling form can be 
checked in Fig. [2] comparing the exact expression (11) to 
the analytical formula (12) for several values of the ratio 
Cg/C of the chain capacitances. Formula (12) predicts 
the rapid decay of correlations Xj = — Z(E 0 /A R ) 2 /j 2 
at large distances. This result is in full agreement with 
the known asymptotics of the magnetic Kondo screening 
cloud in electronic systems. This leads to the identifi¬ 
cation of Eq/A r as the dimensionless Kondo length (or 



FIG. 2. Universal Josephson-Kondo cloud at intermediate 
values of the fine structure constant a = 0.1. Here the ex¬ 
act expression (111 is computed for several choices of the 
capacitance ratio in the chain, C g /C = 10 -4 (dot-dashed), 
Cg/C = 10 -3 (dotted), C g /C = 10” 2 (dashed), C g /C = 10" 1 
(solid), and compared to the analytic scaling form (121. The 
computation is done for Ejd = Eq = 0.01. Vertical lines 
at distances j * = \]C/C g = {VTO, 10, \/T000,100} indicate 
the cross-over between non-universal short-range and univer¬ 
sal long-range behavior. The thin straight line indicates the 
long-distance asymptotic behavior —Z\j = {Eo/A r) 2 /j 2 . 


Kondo temperature A R ). The 1/j 2 behavior may be vis¬ 
ible in chains of several thousand sites, that can be fabri¬ 
cated with modern lithographic techniques. More inter¬ 
estingly however, the decay at intermediate distances, in 
the range yJC/C g < j < Ej/A R , is much slower, attest¬ 
ing to the stronger coupling of the qubit to the Josephson 
elements in its closer neighborhood. It is straightforward 
to see that the existence of these long-range spatial cor¬ 
relations are intimately connected to the strong entan¬ 
glement between the qubit and its environment. Indeed, 
a full polarization of the qubit into the up state, even 
dressed by its cloud of oscillators, leads to short range 
correlations that do not extend beyond the second site 
of the chain. In addition, the long-range correlations of 
the Josephson-Kondo cloud also reveal the non-linearity 
of the Josephson element coupling site i = 0 and site 
* = 1. A purely harmonic chain (see Supplementary Ma¬ 
terial [39]) leads to much faster decaying correlations. 
The Josephson-Kondo cloud is therefore an experimen¬ 
tally accessible hallmark of the strong coupling between 
a two-level system and its macroscopic environment. 

Regime of large fine structure constant. For circuit 
impedances close to the quantum value h/e 2 , the matter- 
light interaction parameter a becomes of order one, and 
the wavefunction ® is in principle not sufficient. It 
was shown recently |37] that the non-trivial many-body 
ground state at large dissipation can be captured by a 
systematic expansion in terms of coherent states, that 
simply generalizes the previous ansatz to a superposi¬ 
tion of M cs coherent states |J/) = Pn |^ n ) , where 
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FIG. 3. Universal Josephson-Kondo cloud at increasing values 
of the dissipation strength a. Fully converged curves (solid 
lines) were obtained with M cs = 7 coherent states. The curves 
were obtained at the following values of the Kondo energy: 
Ar = 4x 1CU 2 I5 o at a = 0.1, A# = 4x 10~ 3 .Eo at a = 0.5 and 
Ar = 5 x 10~ 4 I?o at a = 0.7. In each instance this roughly 
corresponds to Ejd ~ \J ( 2e) 2 Ej/C/20. The M ca = 7 results 
are compared to the M cs = 1 lineshape (121, at the same 
Kondo scale. Lines a ~ 1/j and b ~ 1 /j 2 are there to guide 
the eye. 


each I'I'n,) is a coherent state of the form © with vari- 
ationally determined amplitudes p n and displacements 
4 n) . As M cs is increased, T) rapidly converges to the 
exact ground state. Our technique, based on the ex¬ 
act many-body state, provides a powerful alternative to 
the usual numerical renormalization group calculations 
of the Kondo cloud [36j HSj, which in contrast, requires 
a distinct computation for each spatial distance con¬ 
sidered. The resulting correlation function \j for i n ~ 
creasing values of the dissipation strength a is shown 
in Fig. [3] The cloud polarization \j still decays quadrat- 
ically ~ —Z(E 0 /A R )/j 2 at the largest distances, but 
with a Kondo energy A# that is renormalized down¬ 
wards from the value predicted by (101 (our data at 
a = 0.7 shows a renormalization of A r by a factor 
~ 1/4). At j ~ Eq/Ar, a broad cross-over is visi¬ 
ble. In the intermediate distance universal regime, i.e. 
yJC/Cg <C j <C E 0 /A Rl the correlation function \j de¬ 
cays slower than 1/j. At strong dissipation (large a), our 
results are consistent with the perturbative prediction 
Xj ~ — l/[jln(j) 2 ], but for the realistic system size con¬ 
sidered, it is difficult to achieve the extreme separation 
of scales between sJC/C g and E 0 /A R that is required to 
identify this lineshape unambiguously. We compare the 
converged results with M cs = 7 coherent states to the 
analytical M cs = 1 result (12), but with the true Kondo 
scale Ar, rather than the M cs = 1 estimate (10). It 
is very interesting that, as the dissipation strength in¬ 
creases, a noticeable difference develops at intermediate 
distances between the line shapes predicted by the simple 
M cs = 1 approximation and the fully converged M cs = 7 


result. These differences originate from subtle entangle¬ 
ment within the states in the bath |12j . Indeed, contrary 
to the simple Ansatz © that is based on a single (multi- 
mode) coherent state, the environmental wavefunction in 
the multiple coherent state expansion cannot be factor¬ 
ized. 

Despite the quantitative differences in the intermediate 
distance universal regime, it is very surprising how accu¬ 
rate the simple Ansatz © turns out to be, highlighting 
the unexpected simplicity of the many-body ground state 
of such a complicated and widely investigated Hamilto¬ 
nian as the spin-boson model ([2|. Our results indicate 
that strongly interacting quantum problems display a 
hidden structure that is more generic than the specific 
model studied in this Letter, since fermionic impurity 
problems should show similar properties. With the rapid 
advances in circuit QED architecture, we believe that 
future experiments following our proposal to probe the 
Josephson-Kondo cloud will be able to reveal the sim¬ 
ple correlations that are at the heart of these many-body 
states. 
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by the National Research Foundation of South Africa 
(Grant Number 90657). 


[11 Y. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. Phys. 
73 , 357 (2001). 

[2] R. J. Schoelkopf and S. M. Girvin, Nature 451 , 664 (2008). 

[3] D. Zueco, C. Fernandez-Juez, J. Yago, U. Naether, B. Per- 
opadre, J. J. Garcla-Ripoll, and J. J. Mazo, Supercond. 
Sci. Technol. 26 074006 (2013). 

[4] P. Jung, A. V. Ustinov, S. M. Anlage, Supercond. Sci. 
Technol. 27 , 073001 (2014). 

[5] T. Niemczyk, F. Deppe, H. Huebl, E.P. Menzel, F. Hocke, 
M.J. Schwarz, J.J. Garcla-Ripoll, D. Zueco, T. Hummer, 
E. Solano, A. Marx, and R. Gross, Nat. Phys. 6, 772 
( 2010 ). 

[6] P. Forn-Dlaz, J. Lisenfeld, D. Marcos, J. Garcla-Ripoll, E. 
Solano, C. Harmans, and J. Mooij, Phys. Rev. Lett. 105 , 
237001 (2010). 

[7] O. Astafiev, A.M. Zagoskin, A.A. Abdumalikov, Y.A. 
Pashkin, T. Yamamoto, K. Inomata, Y. Nakamura, and 
J.S. Tsai, Science 327 , 840 (2010). 

[8] A.A. Abdumalikov, O. Astafiev, A.M. Zagoskin, Y.A. 
Pashkin, Y. Nakamura, and J.S. Tsai, Phys. Rev. Lett. 
104 , 193601 (2010). 

[9] E. Sanchez-Burillo, D. Zueco, J. Garcla-Ripoll, L. Martn- 
Moreno, Phys. Rev. Lett. 113 , 263604 (2014). 

[10] K. Le Hur, Phys. Rev. B 85, 140506 (2012). 

[11] M. Goldstein, M. H. Devoret, M. Houzet, and L. I. Glaz- 
man, Phys. Rev. Lett. 110 , 017002 (2013). 

[12] S. Bera, S. Florens, H. U. Baranger, N. Roch, A. Nazir, 
and A. W. Chin, Phys. Rev. B 89, 121108(R) (2014). 

[13] A. C. Hewson, The Kondo Problem to Heavy Fermions 
(Cambridge University Press, Cambridge, UK, 1993). 











5 


[14] A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. 
Fisher, A. Garg, and W. Zwerger, Rev. Mod. Phys. 59, 1 
(1987). 

[15] U. Weiss, Quantum Dissipative Systems (World Scien¬ 
tific, 1993). 

[16] K. Le Hur, Ann. Phys. 323 , 2208 (2008). 

[17] Yu. Bomze, I. Borzenets, H. Mebrahtu, A. Makarovski, 
H.U. Baranger, and G. Finkelstein, Phys. Rev. B 82 , 
161411 (2010). 

[18] D. Y. Baines, T. Meunier, D. Mailly, A. D. Wieck, C. 
Bauerle, L. Saminadayar, P. S. Cornaglia, G. Usaj, C. A. 
Balseiro, and D. Feinberg, Phys. Rev. B 85, 195117 (2012). 

[19] J. E. Gubernatis, J. E. Hirsch, and D. J. Scalapino, Phys. 
Rev. B 35, 8478 (1987). 

[20] V. Barzykin and I. Affleck, Phys. Rev. Lett. 76 , 4959 
(1996). 

[21] I. Affleck, in Perspectives of Mesoscopic Physics: Dedi¬ 
cated to Yoseph Imry’s 70th Birthday, A. Aharony and O. 
Entin-Wohlman (eds.) (World Scientific, Singapore, 2010). 

[22] A. Holzner, I.P. McCulloch, U. Schollwock, J. von Delft, 
and F. Heidrich-Meisner, Phys. Rev. B 80, 205114 (2009). 

[23] C.A. Brisser, G.B. Martins, L. Costa Ribeiro, E. Vernek, 
E.V. Anda, and E. Dagotto, Phys. Rev. B 81, 045111 
( 2010 ). 

[24[ A.K. Mitchell, M. Becker, and R. Bulla, Phys. Rev. B 84, 
115120 (2011). 


[25] I. Affleck and P. Simon, Phys. Rev. Lett. 86, 2854 (2001). 

[26] P. S. Cornaglia and C. A. Balseiro, Phys. Rev. Lett. 90 , 
216801 (2003). 

[27] T. Hand, J. Kroha, and H. Monien, Phys. Rev. Lett. 97 , 
136604 (2006). 

[28] R. G. Pereira, N. Laflorencie, I. Affleck, and B. I. 
Halperin, Phys. Rev. B 77 , 125327 (2008). 

[29] J. Park, S.-S. B. Lee, Y. Oreg, and H.-S. Sim, Phys. Rev. 
Lett. 110 , 246603 (2013). 

[30] J. Bauer, C. Salomon, and E. Dcmler, Phys. Rev. Lett. 
Ill, 215304 (2013). 

[31] I. Affleck, L. Borda, and H. Saleur, Phys. Rev. B 78 , 
139902 (2008). 

[32] T. Heinzel, Mesoscopic Electronics in Solid State Nanos¬ 
tructures (Wiley, Weinheim, Germany, 2007). 

[33] V. J. Emery and A. Luther, Phys. Rev. Lett. 26, 1547 
(1971). 

[34] R. Silbey and R. A. Harris, J. Chem. Phys. 80, 2615 
(1984); 

[35] R. A. Harris and R. Silbey, J. Chem. Phys. 83, 1069 
(1985). 

[36] L. Borda, Phys. Rev. B 75, 041307(R) (2007). 

[37] S. Bera, A. Nazir, A. W. Chin, H. U. Baranger, and S. 
Florens, Phys. Rev. B 90 , 075110 (2014). 

[38] B. Lechtenberg and F. B. Anders, Phys. Rev. B 90 , 
045117 (2014). 

[39] Supplementary Material document. 


Supplementary information for “Josephson-Kondo screening cloud in circuit quantum 

electrodynamics” 


We present here for completeness some intermediate steps of the derivation performed in the main text. 


DERIVATION OF THE EFFECTIVE SPIN-BOSON HAMILTONIAN 
Starting model and definitions 

We start with the Hamiltonian H = Hq + Hi, where the quadratic part Hq describes the linearized semi-infinite 
Josephson junction chain: 


Hn = - 


, +oo 

1 \ ' 


2 ^ 

m,n—0 


Pm 




1 +oo 

Pn + o %mVmn%n ? 

i Z *—' 

m,n—0 


(SI) 


with p m = n m — n m , x m = 4> m and [x m ,p n \ = i5 mn . The impurity located at the end of the chain contains the only 
source of non-linearity of the problem, encoded in the Josephson energy: Hi = —Ej d cos(xq — xi). The matrix 


M = 


1 

J2ey 


( C gd + C d 

-c d 


-c d 

G g + C d + C -C 

-C Cg + 2 c -C 

-C Cg + 2C -C 


\ 


(S2) 




V 
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contains the set of capacitances of the whole system, whereas the matrix 


/ ° 
0 


V = Ej 


0 

1 -1 

-1 2 -1 

-1 2 -1 


\ 


(S3) 


V ••/ 


describes the Josephson couplings in the leads only. The indexing convention is such that the upper-left entries of 
these matrices are labeled by M 00 and Voo- One would like to express H 0 as a bath of independent oscillators that 
are linearly coupled to the impurity at site 0. This is achieved by expressing the set of harmonic variables in terms of 
bosonic operators: 


Pm ~V2 J 

r+oo 

1 dk 0 m k(b k + &£) 

0 

m = 1,2,.. 

., +00 

(S4) 

Xm = v^J 

/*+°° 

d/c £mk(Pk 

0 

TO = 1, 2, . . 

., +0O 

(S5) 


where £ mk and 9 m k are real numbers. Note that we adopted the unusual (but legitimate) convention of ascribing the 
momentum variable (and not the position variable) to the sum of creation and annihilation operators. Since the qubit 
is capacitively coupled to its environment through the momentum po, this choice will result in the standard form of 
the spin boson Hamiltonian. We also stress that the coordinates (xo,po) at the end of the chain are not part of this 
decomposition, which will faciliate the later projection onto the two-level system. A first important relation follows 
from ( S4||S5 ) by imposing the commutation relation [x m ,p n ] = i5 m ^ n and assuming standard algebra of the bosonic 
operators, namely [b k ,b k> ] = 0, [b k ,b\,] = S(k - k')\ 


/» + OO 

I dk ^rnk^nk — 
JO 


(S6) 


Normal mode diagonalization of the chain 


Let’s focus first on the chain part of H 0 which contains the sites m = 1,2,..., +oo only, and impose that it is 
diagonalized by the bosonic operators: 




r+oc, 

dk uj k b 1 ,,b i 


11 

i+fc’ 


(S7) 


m,n^0 


m,n^0 


with ujk the positive frequencies of the normal inodes. Using the decomposition ( S4||S5 ), this provides several 
constraints on the set of unknown parameters: 


m,n^ 0 


m,n^ 0 


@mk 

(*■*) 

@nk' ^mkVmn^nl' 

mn 

dmk 


@nk' H - Cmk^mn^nl' 

mn 


(58) 

(59) 


One can then take the sum and differences of Eq. (S8|S9l, and integrate them respectively with J 0 +oc <ifc £ pk i and with 


fo°°dk 0 p k'. Using the orthogonality condition JS6L one readily gets: 


^ ^ @mk ( M ) — Mk^pk 

z ' V / mv 


m^O 


U CmfeUmp — tdk@pk’ 
m^O 


(510) 

(511) 


In order to write the above equations as full matrix products, we extend the variables to the range m = 0,1,..., +oo, 
and define the vectors = (£ofc, £ifc, &k, ■ ■ •) and @k = ($ofc, ®ik, @ 2 k, ■ ■ •)■ Noting that Vq p = 0 and with the condition 
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9q k = 0, we can write: 

M -1 0fc = (S12) 

v£k=Uk0k- (S13) 

This results in an eigenvalue equation M~ x V^ k = c4Cfc, which can be rewritten in the explicit form (that avoids 
inverting explicitly the matrix M ): 

v£k=oj 2 k M£ k . (S14) 

This gives a set of conditions, for the different cases m = 0, m = 1, and m > 2: 


Co k = 77—r~7rCi k 


Cgd + Cd 


&k = 


1 - 




(2 e) 2 Ej-Cul_ 


Ci fc 


, _ 0 (2e) 2 ^-(C + C' fl /2K, 

Wl - fc (2 e) 2 Ej~CLo 2 k ?mfe 


(515) 

(516) 

(517) 


with C e ff = + CgdCd/{C g d + C'd)■ Now we exploit the fact that the Josephson chain is uniform except at the first 

two sites. We thus introduce scattering states that we parametrize as follows: 


Cmfc = N k cos [k(m - 1/2) -4] m = 1 , 2, ..., +oo 


(S18) 


with amplitude N k and phase shift 5 k to be determined. Inserting the above expression into Eq. (S17), we readily 
find the eigenfrequency: 


2 4(2e) 2 £j sin 2 (fc/2) 

“ C g + 4C sin 2 (fc/2) ' 


The phase shift is then determined from condition (S16), which reads: 

CeS^i 


COS 


3 fc 


- 4 = 


l - 


(2 e) 2 Ej - C4 2 


COS ( - - 4 


After some trigonometric manipulation and using Eq. |Sl9| , one finds the simple condition: 

C„ 


which is solved explicitly as: 


sin (4) = 


4 = arctan 


1 - —sin (4 - fc), 

C c ff 


0 


_ 3l.\ 

C'eff J 


sin(fc) 


.( 1_ &) COs(fc)_1 - 


(S19) 


(S20) 


(S21) 


(S22) 


It is useful to establish an inverse relation between the harmonic variables and bosonic operators. Using (S6), it is 
straight toward to check that the relations (S4]|S5) can be inverted by: 


bfc — r— ^ ^ l^mkPm 'i@mk%m\ 
v m= 1 

(S23) 

1 °° 

z r- ^ ^ [^mkPm H - l i@mk3'm\ • 

^ 771=1 

(S24) 


One is equipped now to determine the normalization factor N k , which follows from the commutation relation [b k , &£,] = 
S(k — k'). This relation is equivalent to the condition 9 k • Cfc' = S(k — fc'), which can be rewritten as: 


/ ' CmfcUmnCnfc' — (i(fc fc )• 

m,n^ 0 


(S25) 






















From the explicit expression of V and the recursion relations ( S16|S17 1, one gets: 

fl _ mi 

9 "£l k€lk' + CmfcC, 


S(k — k') = 4—— sin 2 (fc/2) 

Wfc 


Cn 


mk' 


m—1 


(S26) 


Finally, using expression (S18) and the following standard algebraic identities X) m=1 cos(mfc) = 7r< ^(^) ~~ 1/2 and 
X^m=i s i n (mfc) = l/[2tan(fc/2)], one finds: 


5(k — k') = 2 tt —— sin 2 (fc/2 )N k S(k — k') + Finite Terms, (S27) 

^k 

where we have singled out the delta-function contribution. Let us prove that the finite terms are actually zero, which 
implies the choice of normalization: 


N k = 


Wfc 


sin(fc/2) Y 2-7T Ej 

For this purpose, let us show that 0 k ■ £ k > = 0 for all k ^ k'. The proof relies on equations (S12||S13): 


u 2 J k ■ £ k , = e k ■ M-'vgv = & = Uj 2 k 0 k ■ 


(S28) 


(S29) 


Since ui k is a monotonous function of k, the scalar product 9 k ■ £ k i indeed vanishes for k ^ k'. 


Expression for the full Hamiltonian 


The impurity part of the quadratic Hamiltonian (SI) related to the site m = 0 reads: 

+oo 


Hr = \ 


(m ') pl+Po^(M p m . 
V /OO ' V / Om 


(S30) 


m—1 


Using Eq. (S4) and Eq. (S10), we finally obtain an exact and remarkably compact expression for the full Hamiltonian 
in terms of the normal modes and of the variables at the first two sites of the Josephson chain: 




(M p 2 0 + [dk u> k b\b k + (dk u k t,ok{b k + b\) - E Jd cos(a;o - aq). (S31) 

' 'oo J q v 2 J o 


To get rid of the factors containing aq, we perform a unitary transformation H = UHU', with U = e lp ° Xl . The 
transformed Hamiltonian reads 


1 


H = pi + I dk u k b\b k - p 0 I dk ui k ^ 0 k(bl + b k )-^ [e lx ° +e “°] 


1 




Ejd 


2m 


V2 


where 


1 

m 


(*■') 


oo 


+ 1- 


2 C d 


Cd + C d g J J 0 


dk u k (l k . 


(S32) 


(S33) 


Now we assume that we are at a degeneracy point of the uncoupled impurity, and project the impurity part of the 
Hilbert space onto the two degenerate states. This results in p 0 —t er z /2 and ( e lx ° + e~ lx °) —>• a x . It produces a 
spin-boson Hamiltonian (where again, we drop constant terms) 


H ~ / ui k b\b k - 


For k <C sJCg/C , 


i c dg r 

k x/ 2 C dg + C d Jo 


co k ~ E 0 k 
Ci fc — 2 


dkw k £ lk (bl +b k ) a C - Ej d °C- 


Eo 


2nEjk 


(534) 

(535) 

(536) 
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with E 0 = \J (2 e) 2 Ej/C g . We can thus identify the Ohmic spin-boson model parameter 


C, 


dg 


1 


27T \ Cdg + Cd 


Eo 

E~j' 


(S37) 


From the point of view of the impurity, a is the effective fine structure constant when the photonic modes of the 
vacuum are replaced by the plasmonic modes of the superconducting environment. This interpretation is confirmed 
by noting that the fine structure constant in vacuum a V ac. = (e 2 /2 h)Z va , c , (with e the electron charge and h Planck’s 
constant) is fixed to the small 1/137 value by the vacuum impedance Z vac . = y/Jiojco — 377 Q. In transmission 
lines, an effective fine structure constant can thus be defined as a = (e 2 /2 h)Z from the environmental impedance 
Z, typically related to a ratio from line inductances L and shunt capacitances C , Z = yjL/C . In superconducting 
circuits, the Josephson inductance is given by L = (h/ 47re) 2 l/ Ej , so that we indeed recover Eq. (S37) up to numerical 
and geometrical factors. 


Connection between the magnetic Josephson-Kondo screening cloud and the Kondo screening cloud 

Here we explain the mathematical correspondence between the Kondo screening cloud and the Josephson-Kondo 
screening cloud. For a spin 1 /2 magnetic moment at x = 0, quenched in a one dimensional electron gas, the longitudinal 
Kondo screening cloud is defined as 


xii( S ) = 4<5r p 5: i (^)>. 


(S38) 


Here Sl mp is the ^-component of the impurity spin, and S^(x) is the ^-component of the electron spin density at x. 
The cloud can be decomposed into 


X^{x) = Xq(:e) + cos(2k,Fx)x\ kF (x), (S39) 

where Xq(x) and X^ kF (x) vary slowly on the scale of the Fermi wavelength 2’K/kp. Using the well-known mapping 
between the Kondo model and the spin-boson model, PE] the component (x) can be expressed as 


aXl (x) = 



where 


X q {x) = y/qe q/2 cos (^j . 


(S40) 


(S41) 


Here 1/a is the ultra-violet scale in the problem (of the order of /cp), q is dimensionless, and the expectation value is 
with respect to the ground state of the spin-boson model. Comparing this expression to the expression (7) for i n 
the main text, we see that the two integrands have the same low q behavior. Thus, the 0 ItF component of the Kondo 
cloud for x a, corresponds to the Josephson-Kondo cloud for j yJCjC g . Explicitly, the correspondence is 


X ^j~l+Sl (S42) 

a 2 


At distances where the correspondence is accurate, the profile of the cloud is universal. This can be seen analytically 
in the regime of small to moderate dissipation, where the exact Josephson-Kondo cloud is given by Eq. (11) in the main 
text. As pointed out there, it depends on all the microscopic parameters of the device. However, for j sufficiently larger 
than y/C/Cg, but not necessarily larger than 1/A#, the integrand of Eq. (11) in the main text can be approximated 
as 


cos[fc(j — 1/2 + SI)] 


E 0 k + A# 


(S43) 
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due to rapid oscillations when k > l/j. For the same reason, the upper bound of the integral can be extended from 
7r to oo. Under this approximation, the integral evaluates to 


Xj = -ZRe { e" i ^ ( - 7 ' _1/2+,50 r 


-1/2 + 61) 
■C'O 


(S44) 


where r(a, z) is the incomplete Gamma function and Z is as defined in Eq. (13) in the main text. We have said that 
under the assumptions we make about capacitances, SI is small. In the same regime, the Kondo scale is small 
relative to the ultra-violet scale E 0 . This implies that \j varies slowly on the scale of Si — 1/2, and we can therefore 
replace j — 1/2 + SI —> j in the above equation. In this way, we obtain Eq. (12) in the main text. 


The harmonic chain 


In the main text we argued that the non-trivial physics we investigate is generated by the anharmonic term 
—Ejd cos(<po — <f>i) in the Hamiltonian (1). Specifically, we want to show that without the anharmonicity, correlations 
decay far more rapidly than with the anharmonicity. In order to justify these statements and make them more precise, 
we study here the harmonic chain. That is, we consider the regime E Jd (2 e) 2 /(C d + C gd ), where the cosine term 
in the Hamiltonian (1) can be expanded to quadratic order in the phase difference </>o — 0i, yielding a harmonic chain 

Hh = ( Ui ~ (^ _1 ) .. ( n l - n i) + - <t>i+ 1) 2 + ^ ~ M 2 - ( S45 ) 

i,j -0 lJ i=l 


The Hamiltonian is again diagonalized by dehning bosonic operators 4 and b\. , for k € [0,7r], such that 


4 


1 

71 


7^. [£jk( n j n j) 
1=0 


(S46) 


As before, 4fc is obtained from the classical equations of motion. The presence of the quadratic Josephson coupling 
between sites i = 0 and i = 1 modifies the expressions for as follows. The relation between £o k and £ik now reads 


^0 k (1 Ek^^iki Ek 


ulC gd 

L 0 2 k (C gd + C d )-(2e) 2 E Jd 


(S47) 


For j > 1, the form of £jk is still 


tjk = N k cos [k(j - 1/2) - 4], 


(S48) 


In the last two equations, the from of u>k and Nk are unchanged from what they were before, but the phase shift 4 
is obviously affected by the Josephson coupling between sites i = 0 and i = 1. It is now given by 


tan(4) = 


sin k 


cos k — 1 — 


C gi (l-F k ) 


(S49) 


The relationship between 4 and 4 is still given by (S12) and (S13), but in (S13), one has to remember to include 
the extra Ej d couplings, i.e. Vqo = Ej d , V\\ = Ej d + Ej, and Vqi = Vio = —Ej d . The ground state of the harmonic 


chain (S45) is the bosonic vacuum |0) defined by the condition that 4 |0) = 0 for all k in [0, n\. 
For the correlator \j °ne then obtains 

Xj = (0| (n 0 - n 0 ){n 3 - %) |0) 

1 r 

= 2 J ^0k@jk 

e r f 

= -- / dk — cos(fe/2 — 4) cos[k(j — 1/2) — 4]> 

7r Jo 


(S50) 


the last line being valid for j > 2. Note that, while F k diverges at the point k = fc* such that to 2 = (2 e) 2 Ej d /{C gd + 
C d ), the product F k cos(fc/2 — 4) remains finite at k = so that the integrand in the last line of (S50) is a smooth 
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j 


FIG. SI. The correlation function \j for the harmonic chain, calculated from Formula (S501. The dashed curve shows the result 
for weak dissipation (large C), and was calculated for Cd = C g d = C g , Ejd = 10(2 e) 2 /C g , and Ej = 10(2e) 2 /(2C + C g ). For 
the anharmonic impurity system (i.e. a much smaller value of Ejd), the remaining parameter values would have corresponded 
to a = 0.035. The solid curve shows the result for strong dissipation (small C), and was calculated for C = Cd = C g / 100, 
C g d = C g , and Ejd = Ej = 10(2e) 2 /C g . For the anharmonic impurity system (i.e. a much smaller value of Ejd), in this case, 
the remaining parameter values would have corresponde d to a = 0.5. For both the solid and the dashed curve, a thin solid line 
indicates 1 / j 2 with the prefactor predicted by Formula (S52). 


function of k. At sufficiently large j, the slowly varying factor FkCOs(k/2 — 5k)/uk can be evaluated to linear order 
in k, yielding 

Xj - ~ -J7 E~yT f dkk co #0' - V 2 ) - 4], 
n (2e) Jo 

= - E ,° C , 9 2 d j [ dk sin [k(j ~ 1/2) - 4], (S51) 

^ 44 Jo 

At large j, the phase shift 4 dephases the contribution from evaluating the integral without the phase shift at k = n. 
Replacing j — 1/2 —»• j, which is allowed within the accuracy to which we evaluated the integral, we then find 


1 E 0 C g d „ 


— - cos (jk) 

- J Jo 


EpC gd 1 
7r (2e) 2 j 2 ' 


For the anharmonic system that we studied in the main text, we found on the other hand 


(S52) 


v M 2 1 X ^Cgdl_ 

\A R ) f 2A R (C d + C gd ) n(2e) 2 j 2 ' 


(S53) 


At strong dissipation, this results in a correlation function that is larger than in the harmonic case by a factor 
proportional to the inverse of the (exponentially small) Kondo energy. 

Next we ask at what distance the 1 /j 2 behavior predicted in Formula (S52) sets in. For the harmonic chain, the 1/j 2 
decay breaks down below a distance j ~ l/fco, where ko is the value of k where the fc 3 term in the Taylor expansion 
of Fk cos(fc/2 — 8k)/oJk in the integrand of ( |S50 ) becomes important. Under the assumptions we make about the sizes 
of capacitances, it follows that k 0 ~ min {JJUgJC, 1}. Hence, for the harmonic chain, \j ~ 1/j 2 behavior is seen for 
j max{ C/C g , 1}. In contrast, for the anharmonic impurity, we saw in the main text that there is an intermediate 
regime \JC/C g < j <C Eq/Ar, where correlations decay far more slowly. 

We also note the following. In the anharmonic impurity problem, one increases the dissipation strength a to 
values of order 1 by decreasing C, while keeping Ej{2C + C g )/{ 2e) 2 large. This has the effect of decreasing Ar, 
and hence the region of the onset of 1/j 2 decay is shifted to larger j. In the harmonic problem, increasing C tends 
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to have the opposite effect, namely, it shifts the onset of 1 /j 1 2 decay to smaller j. This is illustrated in Fig. SI 
where the exact formula ( |S50 1 and the large j asymptotic formula (S52) for \j i n the harmonic system are compared, 
for fixed Ej(2C + C g )/(2e) 2 , but different C. These results justify the statement we make in the main text that 
the Josephson-Kondo cloud is a hallmark of the strong coupling between a two-level system and its macroscopic 
environment. 
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